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Abstract. BAT99-129 in the LMC is one among a handful of extra-galactic eclipsing Wolf-Rayet binaries known. We present 
blue, medium-resolution, phase-dependent NTT-EMMI spectra of this system that allow us to separate the spectra of the two 
components of the binary and to obtain a reliable orbital solution for both stars. We assign an 05V spectral type to the com- 
panion, and WN3(h)a to the Wolf-Rayet component. We discuss the spectroscopic characteristics of the system: luminosity 
ratio, radii, rotation velocities. We find a possible oversynchronous rotation velocity for the O star. Surprisingly, the extracted 
Wolf-Rayet spectrum clearly shows the presence of blueshifted absorption lines, similar to what has been found in all single 
hot WN stars in the SMC and some in the LMC. We also discuss the presence of such intrinsic lines in the context of hydrogen 
in SMC and LMC Wolf-Rayet stars, WR+O binary evolution and GRB progenitors. Altogether, BAT99 129 is the extragalactic 
counterpart of the well-known Galactic WR binary V444 Cygni. 



Key words. LMC - massive star - Wolf-Rayet star - O star - Wolf-Rayet: BAT99-129 - 
classification - spectroscopic binary - eclipsing binary - Gamma-ray burst progenitor 
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1. Introduction 

The Wolf-Rayet (WR) st ar BAT99-129 (a.k.a Brey 97, 
see iBrevsacher et alJ [1999, hereafter BAT129) has been 
discovered to be a sho rt -period eclipsing binary by 
iFoellmi. Moffat. & Guerrerol J2003b). According to these 
authors, it consists of a nitrogen-rich WN4 Wolf-Rayet compo- 
nent and a companion of unknown spectral-type. Along with 
BAT99-19 (Brey 16) in the LMC and HD 5980 in the SMC, 
BAT129 belongs to the very small group of extra-galactic 
eclipsing Wolf-Rayet binaries. 

Absorption lines are visible in the spectrum of BAT129, es- 
pecially in the blue part of the optical range (<4000A), while 
WR emission lines strongly dominate everywhere else. The ori- 
gin of these absorption lines is unclear: they are either pro- 
duced by an O-type companion, and/or they originate in a 
WR wind, as has been found in all the single early-type WN 
(WNE) stars in the SMC and, to a lesser extent, in the LMC 
jFoellmi. Moffat. & Guerrerd2003albllFoellrrul2004l) . 

In the SMC about half of the WR population consists of 
hot, single, hydrogen-containing single WN3 and WN4 stars 
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(classified WN3-4ha or (h)a). For such hot spectral types, sig- 
nificant amounts of hydrogen were not expected, since such 
hot WN stars are believed to have peeled off their outer Fi- 
nch layers by a very strong and optically thick wind, combined 
with internal convective mixing to expose core-processed ma- 
terial. Contrary to expectations, blueshifted absorption lines of 
HI and Hell were detected in the spectra of all single SMC 
WNE stars, and Foellmi et al. argued that they originate in a 
WR wind. Such stars were also found in the LMC, and more 
recently, in our Galaxy, in a region where the ambient metallic - 
ity may resemble that o f the LMC: WR3, classified as WN3ha 
( Ma rchenko et alJl2004l) . Finally, two other Galactic, presum- 
ably single WNE stars clearly show hy drogen in their spec- 
tra: WR128 (WN4) and W R152 (WN3) JCrowther etalll995t 
iNugis & Niedzielskll995T) . 

This paper examines a binary star. Hydrogen is already 
known to occur in the WR component of som e binaries, e.g. 
in the 1.64-day eclipsing system CQ Cep jMarchenko et alJ 
Howe ver, th e auth ors classify the WR component as 
WN6, while ISmithl J 1968b gives a spectral type WN7. These 
two spectral types seem to mostly describe very massive stars 
(see for instance the most massive star known: WR20a in our 
Galaxy, with two stars both of type WN6ha - Bon anos et alJ 
120041 tRauw et alJ 120051) in a core hydrogen-burning transi- 
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tion stage between the O M ain Sequence an d the core helium- 
burning WR stage (see also lFoellmi et all2003bl for a detailed 
discussion of this important point). The question that remains 
to be addressed is: do we detect the presence of hydrogen in the 
early-type WN component of a short-period binary, and if yes, 
could such a presence be explained in the same way? 

Previously, none of our spectra of BAT 129 were of suffi- 
cient quality to detect absorption lines belonging to the com- 
panion and measure its radial-velocities (RVs), thus obtaining 
a reliable orbital solution for both components and helping to 
separate the two spectra, as well as to resolve the question of 
hydrogen intrinsic to the WN component. 

We present in Section 2 new medium-resolution NTT- 
EMMI spectra of BAT129. We discuss in detail in Section 3 
the iterative procedure we used to separate the spectra and find 
a reliable orbital solution. In Section 4, we assess the presence 
of hydrogen in the WR spectrum and examine all the results on 
this binary star: spectral types of both components, luminosity 
ratio, radii, mass-loss rates, masses. We discuss in Section 5 
the relevance of the presence of hydrogen to the WR classifica- 
tion, WR+O binary evolution and GRB progenitors. We draw 
conclusions in Section 6. 



2. Observations and Data 

Our spectroscopic dataset of BAT129 is composed of three 
parts. The first part (binary discovery) was published by 
iFoellmi et al.l d2003bl) . to which the reader is referred for a de- 
tailed description. There one can also find the MACHO data 
used to obtain the light curve of BAT 129. The second part con- 
sists of 21 medium resolution spectra spread over a few days, 
obtained with the ESO Multi-Mode Instrument (EMMI) at the 
New Technology Telescope (NTT) on La Silla, Chile. The third 
part is composed of 12 spectra with identical characteristics to 
the second dataset, but taken six months later. Tablensumma- 
rizes the characteristics of each dataset. 

For the new EMMI data, the Blue Medium-Resolution 
Spectroscopy mode (BLMD) was used, with Grating#3 at a 
central wavelength of 4150 A, providing a wavelength cov- 
erage of 3920-4380 A and a dispersion of 0.45 A/pixel. The 
l"-wide slit provided a resolving power of 3500 (resolution of 
2.6 pixels). The typical 30-min exposure resulted in S/N ratio 
between 50 and 90. The spectra were reduced using standard 
IRAF tasks and shifted to the heliocentric rest frame. Using 
continuum regions common to all EMMI spectra, the contin- 
uum was normalized to unity. The mean spectra of dataset 1, 
and datasets 2 and 3 together, are shown in Fig. [2 It can be 
seen that the new spectra allow us to clearly resolve the emis- 
sion and absorption components. 



3. Analysis 

We performed an iterative analysis on the EMMI spectra, that 
leads to a new and reliable orbital solution as well as separa- 
tion of the individual spectra of the components. The general 
appearance of the extracted spectra, and therefore the shapes 
and strengths of features found in them, are very sensitive to 
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Fig. 1. Top panel: Mean spectrum of data from Foellmi et al. 
2003b. Bottom panel: EMMI mean spectrum. In this spectrum, 
the emission-line and absorption-line profiles are clearly re- 
solved. 



the process of separation. We describe below in detail the anal- 
ysis and the decisions we took in order to reach a reliable result. 

3.1. Photometric period, ephemeris and continuum 
correction 

Before computing the orbital phases for the EMMI spectra, 
we recalculated the ephe meris (p eriod and zero point Eq) pub- 
lished by Foellmi et al. (2003b). The published values were 
based on the orbital solution obtained for RVs of emission 
lines only. In our new approach, we tried to minimize any 
influence of wind-wind collision (WWC) effects on the RVs. 
Furthermore, emission-line RVs of WR stars may not be a good 
indicator of the systemic RV, although orbital motion may not 
be a problem per se. This point is also discussed below. 

A period search in the combined R and B MACHO data 
led us to a period of 2.7689 + 0.0002 days, which is equiva- 
lent to the published value. We note that with the I-band OGLE 
data (ID: LMC.SC18 122q46/OGLE054148.68-7035 31.0; see 
also lFoellmi et alJl2003bl) . IWvrzvkowski et all J2003) obtain a 
value of 2.76895 days. To determine the zero point (i.e. time of 
primary minimum) Eq, we applied the bisector method on the 
deeper eclipse light-curve minimum, and obtained the follow- 
ing value: E [HJD] = 2448843.8935 + 0.0007. The resulting 
light curve folded in phase is shown in Fig. [2 

In order to be used in the separation procedure, the EMMI 
spectra must be corrected for the variable continuum level. 
We smoothed the lightcurve by taking mean values from 0.01 
phase bins, converted magnitudes into fluxes and interpolated 
the smoothed lightcurve at the appropriate phases of the RV 
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Dates 


Observatory 


Instrument 


No. of spec. 


Wavel. coverage 


Res. power (2.6 pix resol.) 


1998-2002 


See (1) 


See(l) 


24 


3800-6800A 


-900 


2004, Sep. 22-25, Oct 1st 


NTT 


EMMI 


21 


3920-4380A 


3500 


2005, Mar. 24-26 


NTT 


EMMI 


12 


3920-4380A 


3500 



1 4.7 fh — 1 — 1 — 1 — i — 1 — 1 — 1 — I — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — r= 
14.8 - 




15.2 -". 



0.2 0.4 0.6 0.8 1 

Phase 

Fig. 2. MACHO ligh t curve of BAT129, as described in 
iFoellmi et all J2003bl) . with the new value of E . MACHO 
instrumental magnitudes were transformed to obta in a mean 
magn itude outside the eclipse of 14.92 JBrevsacher et alJ 
1999). The phases at which EMMI spectra were obtained are 
indicated by vertical lines. Full lines indicate spectra from the 
second dataset, whereas dashed lines indicate spectra from the 
third dataset. 



data. Then we computed the ratio between the flux at a given 
phase and the mean flux taken outside the eclipses, and appro- 
priately normalized each spectrum by this ratio. The correction 
was at most 15%. In that manner, the shapes (and equivalent 
widths) of the spectral lines do not depend on the variable con- 
tinuum level. 

3.2. The "old" radial-velocities 

We also decided to re-measure the RVs of the first dataset, 
clearly affected by the traveling WWC emission details. New 
measurements were made using the bisector method on the 
strongest and best defined line: Hell ,14686. 1 The bottom level 
for the bisectors was fixed at 1.125 (the continuum being at 
unity), in order to avoid blending problems. To minimize the 
influence of the emission component arising in the WWC zone, 
we fixed the top level at 1.5, thus using the lower-third portion 

1 We also assessed the Hell /1541 1 line, but discarded the measure- 
ments because of the relatively low S/N ratio in ~ 50% of the spectra. 



of the line. This proved to be quite reliable, since the results did 
not change significantly if we slightly modified these limits. 

3.3. The preliminary radial-velocities of the absorption 
profiles 

Before measuring RVs in the EMMI spectra, we used the two 
narrow interstellar Calcium lines H & K at 3933 and 3968 A, 
clearly visible in the EMMI spectra, to self-calibrate the RVs 
by co-aligning the IS features in the individual spectra, in RV 
space. The correction was at most 30 km/s for a given spectrum. 
We then removed these lines using gaussian fitting, in order to 
avoid problems during the separation process (see below). 

Simultaneous measurement of the emission and absorption 
profiles in the original EMMI spectra proved to be problem- 
atic. The uncertainties are introduced both by (a) the large- 
amplitude relative motion of the absorption and emission pro- 
files, and (b) the traveling emission component produced in the 
WWC zone. We therefore chose at that point to measure the 
absorption profiles only, since they are relatively well-defined 
in each individual spectrum. 

We measured positions of the absorption profiles via least- 
square fitting of Gaussian profiles (the IRAF task "splot") for 
the following lines: Hell ,14025, U6, Hell ,14200 and Hy. A 
weight of 0.5 was assigned to the relatively weaker absorptions 
at ,14025 and ,14200, compared to 1 .0 for the rest. We thus com- 
bined these four lines into a single set of absorption-line RVs. 2 

3.4. The first-iteration orbital solution 

Combination of the newly obtained RVs of the emission lines 
(dataset 1) with the RVs of the absorption lines in EMMI 
spectra leads to a complete set of spectroscopic orbital pa- 
rameters. We applied the algorithm described in detail by 
Bertiau & GrobbenJ i 19691) to obtain the solution. We fixed the 
period to the photometric value. Since the derived eccentricity 

2 The situation is slightly more complicated since Hell /14025 is in 
fact a blend of the lines Hell ,14025 .60 and Hel ,14026. 191. Moreover, 
the Balmer lines H<5 and Hy are blended with Hell lines as well. In 
order to be perfectly precise, the positions of the profiles measured 
on the spectra should be compared to the rest wavelength of the cor- 
responding lines. But in the case of blends, the contribution of each 
individual line is not known. We make here the reasonable hypoth- 
esis that the relative intensity contribution of one individual line to 
a given blend is not changing with time. Therefore, the comparison 
of a profile position to one rest wavelength is relevant and coherent. 
Thus the choice of one rest wavelength to another will only produce 
an offset in the absolute value of the RVs. At that point of our analysis, 
only relative RVs are meaningful. In general, except when specifically 
stated, the rest wavelength that is used is the one of the line specified 
explicitly in the text. 
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Table 2. First iteration of orbital parameters of BAT 129. The 
parameters published by Foellmi et al. 2003b are also shown 
for comparison. 



Paramptpr 

1 £11 dlUt LCI 


Fnellmi et al ?003h 

1 \JC111111 CI £11. Z.UVJU 


± iiifl aiuuv 




7 7687 + 0007 


7 7689 ("fixed") 


K„,„ (Vm 
l^-WR (Kills ) 


+ 7 

J70 X / 


J 1 O HI / 


K (km s -1 ) 




174 ± 7 


e 


0.082 ± 0.022 


0. (fixed) 


E (HJD-2 451 100.0) 


820.47 ± 0.09 


846.4468 ± 0.009 


y (km s" 1 ) 


245 ± 15 


270 ± 5 


wf) 


342 ± 12 


270 (fixed) 


a. WR sin i (R B ) 




17.5+0.5 


a sin i (R 3 ) 




9.5 ± 0.6 


o-tO-CWkms- 1 ) 




39 


o-(0-C) (kms- 1 ) 




17 



did not significantly differ from zero, we subsequently adopted 
e = 0.0 (but see below). 

On the other hand, we allowed the emission-line RVs to ex- 
hibit a phase shift. Although our bisector measurements aimed 
to avoid (at least, substantially minimize) any WWC influence, 
it is still possible that these RVs may show residual WWC or 
other effects, therefore misrepresent ing the genuine orb ital mo- 
tion (see the detailed discussion in Lewis et al.lll993l We al- 
lowed the phase shift to be an additional free parameter, in- 
troducing it in discrete steps of 0.001 days, from -0.2 to 0.2 
days. We thus computed a series of orbital solutions, artifi- 
cially shifting the time of emission-line RVs with these step 
values, and selecting the solution with the minimal rms devia- 
tion between the observed and calculated RVs. Table |2]shows 
the com parison between thi s solution and that previously pub- 
lished dFoellmi et all2003bl) . 

Clearly, the newly obtained RVs provide much smaller Kwr 
amplitude, by about 80 km s _1 . This is due to the newer, better 
emission-line RVs from dataset 1, and resulting higher-quality 
fit of our RVs. The phase shift between the two sets of RVs 
is equal to 0.012 days, i.e. 17 minutes, which is negligible. 
Moreover, the maximum RV difference between this solution 
and a solution with a null phase shift does not exceed 5 km s . 
The two sets of RVs can therefore be considered as being well 
in phase with each other. 

3.5. The first attempt at spectral separation 

We used the calculated RVs in an attempt to separate the 
spectral componen ts in the EMMI data following th e iterative 
method from Dem ers et al.l d2002l) (See also March enko et al.l 
19981). This method consists of the following steps, starting 
with the emission-line component: 

1 . Shift the spectra (always in RV space) to the common rest 
frame of the WR star, using the (fitted) orbital RVs of emis- 
sion lines, 

2. Combine the shifted spectra (possibly weighted by the indi- 
vidual S/N ratios), thus providing a "super WR spectrum", 
where the O-star features have been smeared out. 



3. Shift the resulting mean spectrum back to the original po- 
sitions, one by one. 

4. Subtract the mean spectrum from the original spectra, thus 
obtaining a set of first-iteration, individual spectra of the O 
component. 

5. Using the RVs of absorption lines, shift these first-iteration 
absorption-line spectra to the common rest frame of the O 
star. 

6. Combine them (possibly weighted by individual S/N ra- 
tios), thus providing a mean first-iteration absorption-line 
spectrum. 

7. Shift the mean absorption-line spectrum back to the origi- 
nal positions, one by one (either using the previously cal- 
culated RVs or estimating them via cross-correlation). 

8. Subtract this mean absorption spectrum from each of 
the original spectra, thus providing a better guess at the 
emission-line spectrum and finishing the iteration. 

The process is repeated until a stable solution is found, usu- 
ally after 3 to 5 iterations. The choices to be made, as well as 
sources of uncertainty in the separation procedure are the fol- 
lowing: 

- The weights during the combination of spectra. 

- The possibility of editing by hand the combined spectra. 

- The possibility to compute (via Cross-Correlation = CC) 
new RVs in the currently extracted spectra, in each itera- 
tion, before shifting. 

- The wavelength range used in the CC. 

- The reference spectrum used in the CC. 

- The spectra to be combined into a mean spectrum (all 
phases, undiscretionately or some only at specific phases?). 

After many tests, we converged to the following: 1. We de- 
cided not to use weights when combining the spectra in order 
not to favor any orbital phase, as we are dealing with a rel- 
atively small number of spectra. Most of our EMMI spectra 
have about the same S/N ratio (between 50 and 90). This op- 
tion proved to be of negligible influence. 2. Editing by hand 
the combined spectra proved to be important in order to force 
the process to converge rapidly. The editing was done in order 
to minimize the adverse alteration of the spectrum, by simply 
erasing absorption profiles from the WR spectrum and emis- 
sion details (presumably related to the WWC effects) from the 
absorption spectrum. Tests reveal that this is absolutely neces- 
sary for adequate convergence. 3. When using the CC for com- 
puting new RVs before shifting, we chose the wavelength range 
4000-4370A, i.e. the spectral range comprising the 4 most im- 
portant absorption lines: Hell /14025, H<5, Hell i4200 and Hy. 

4. The reference spectrum used in the CC computations was 
chosen to be the average spectrum from the previous iteration. 

5. Finally, the spectra around phase 0.9 show very clear WWC 
emission line excess, and their profiles considerably influence 
the resulting mean spectrum. They are certainly not representa- 
tive of an unperturbed WR spectrum, and their influence must 
be tested. 

In order to test our separation process, we performed 4 
different separations, called hereafter "loops". In the first two 
loops we used the fitted RVs from the orbital solution above 
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during the first 2 iterations, then switching to the RVs com- 
puted with CC in 3 additional iterations. In the first loop we did 
not take the spectra at phase 0.9 when compiling the averages. 
Note that the 0.9-phase spectra were kept in the process, and we 
obtained extracted components at those phases, but they were 
not used during the combination of the spectra (steps 2 and 6). 
In the second loop we used all the spectra. 

The two other loops used the RVs from the orbital solution 
only, during 5 iterations. No new RVs were computed with the 
CC method, thus final results being a 'black-box' output of the 
iterative process, based on preliminary RVs of the first orbital 
solution. Again, we removed the phase 0.9 spectra when cal- 
culating the mean spectrum in the third loop, and treated all 
phases equally in the fourth loop. 

The resulting spectra of the 4 loops were found to be very 
similar, at first sight. However, significant differences were 
found in the RVs that were computed from them. In fact, since 
we don't know a priori the correct shape of the spectrum of the 
WR and the O components, we have to find a way to discrim- 
inate between these 4 approaches. We used orbital solutions to 
achieve this. 

3.6. Exploring the parameter space of the orbital 
solutions 

As stated in the first orbital solution above, the emission-line 
RVs from dataset 1, and the absorption-line RVs from EMMI 
spectra seem to be well in phase. To check this point again 
with the new absorption-line RVs on extracted O-star spectra, 
for each of the 4 different loops, we computed the RVs for the 
O spectra via CC in the wavelength range 4000-4370 A. Then, 
we looked for an orbital solution taking only these RVs and 
the emission-line RVs of dataset- 1, allowing for a phase shift 
(i.e., again taking the solution with the smallest cr(0 - C)) as 
described above. We found 4 solutions in quite good agree- 
ment with the first orbital solution (see section l3~4l l. with phase 
shifts ranging from -0.067 days (-1.6 hour, 2% of the period) 
to -0.012 days (0.3 hours, 0.4% of the period). This can be re- 
garded as negligible, meaning that these EMMI absorption-line 
RVs obtained on extracted spectra are well in phase with the 
RVs of dataset- 1, and give a stable overall result. These two 
sets of RVs can therefore be combined into a single set, and 
give us confidence on the results of the separation algorithm. 

Second, we looked for orbital solutions taking the three 
sets of RVs, allowing for a phase shift of the emission-line 
RVs from the extracted WR spectra only (since the previous 2 
datasets were found well in phase). Moreover, we assume here 
that the number of spectra we have and the coverage in phase 
they provide are good enough to find a reliable solution, with- 
out additional constraints such as a circular orbit. We assume 
that the best solution is the one which gives the minimal scatter 
between fitted and observed RVs, and the smallest value of the 
eccentricity and the smallest phase shift. For each of the 4 dif- 
ferent loops, we looked for an orbital solution with a free phase 
shift between -0.2 and 0.2 day in steps of 0.001 day and with 
eccentricity as a free parameter. 



For each of the 4 loops, we have two "best solutions", one 
with the smallest scatter cr(0 - C) and one with the smallest 
eccentricity. It appears clearly in our results that the third loop 
(i.e. blind separation of 5 iterations without combining the 0.9- 
phase spectra) provided both smallest values for cr(0 - C) (of 
36.6 and 39.2 kms _I , respectively), while all other solutions 
have a cr(0 - C) above 40 km s _1 . The set of spectra providing 
the orbital solutions of the best quality is therefore that of the 
third loop. 



600 - 



400 



200 



-200 




0.4 0.6 
Phase 



Fig. 3. Final (adopted) circular orbital solution including all 
3 sets of RVs. Squares: RVs of emission lines (dataset 1). 
Triangles: RVs of emission lines (datasets 2 and 3). Circles: 
RVs of absorption lines (datasets 2 and 3). 



Table 3. Adopted circular orbit of BAT129. 



Parameter 


Circular orbit (adopted) 


P (days) 


2.7689 (fixed) 


K WR (kms- 1 ) 


316 ±5 


K (kms- 1 ) 


193 ±6 


c 


0. (fixed) 


E 


2451946.4274 +0.0054 


y (kms -1 ) 


265 ±5 


a WR sin i (R e ) 


17.3 ± 0.4 


a sin i (i? Q ) 


10.6 ± 0.5 


otO-CWkms- 1 ) 


34.0 


a-(0-C) (kms- 1 ) 


13.4 


M WR sin 3 i (M ) 


14 ± 2 


M sin 3 i (M Q ) 


23 ± 2 



It is not too surprising that a "blind" separation loop pro- 
duces the best results. Although 2 blind iterations (in loops 1 
and 2) allow the process to converge rapidly to mean WR and 
O profiles, the next iterations with RVs computed from them 
are still very sensitive to any residuals. 
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When comparing the two "best" solutions of the third loop, 
it is clear that they are very similar, and can be considered as 
equivalent. All orbital parameters are identical within the er- 
rors, except for the phase shift. The solution with the smallest 
cr has a phase shift of 0.127 days (3.05 hours), while the one 
with the smallest eccentricity has a phase shift of 0.076 days 
(1.8 hours). Given the fact that both solutions have a resulting 
eccentricity consistent with an almost circular orbit: 0.022 and 
0.021 + 0.015, respectively, we decided to choose the solution 
with the smallest cr (i.e. the one with the 0. 127-day phase shift), 
and enforce a circular orbit for simplicity. Finally we adopted 
the circular orbit, summarized in Table[3] and shown in Fig. [3] 



gions and then with a minimal number of samples. The further 
steps were performed identically and led to the same sequence 
of decisions between the different orbital solutions. Finally, we 
chose the continuum rectification that provides the best (i.e. 
flattest) continuum on the final, extracted spectra, i.e. the one 
with the maximum continuum samples. 

4. Results 

We discuss here the results of the separation loops, i.e. the re- 
sulting separated spectra, their properties, and the orbital solu- 
tion. 



3.7. The final separation 



4.1. The Wolf-Rayet absorption features 



Using the fitted RVs from the last orbital solution, we per- 
formed an ultimate 'blind' (no CC measurements) separation 
with 5 iterations, with a limited set of spectra to combine (i.e., 
excluding phase 0.9) and no weights applied during the aver- 
aging. The extracted spectra are shown in Fig.0] 



1.1 



0.9 



01 





Extracted spectrum 



3900 



4000 



4100 4200 
Wavelength [A] 



4300 



4400 



Fig. 4. Individual spectra of the WR and O components. The 
WR spectrum has been shifted vertically by 0. 1, while keeping 
the same scale. 

It must be noted here that the shapes of the wings of ab- 
sorption and emission lines are very sensitive to the adopted 
levels of local continuum. Given the broadness of the emission 
lines and their substantial shifts due to orbital motion, it is dif- 
ficult to define adequately wide common line-free regions to be 
used during continuum rectification. The most apparent conse- 
quence of a non-optimal continuum fit is the presence of low- 
intensity emission excesses in the wings of absorption lines in 
the O-star spectrum (see Fig.0}. To alleviate this problem, we 
repeated the whole separation procedure (i.e. from section l3~2l 
to !3.7l twice, changing the approach to the continuum rectifica- 
tion. The continuum fits were done with the same polynomial, 
but once with a maximal number of available continuum re- 



One important detail can be seen in the separated WR spec- 
tra: absorption profiles on the blue sides of the emission lines. 
These features were always present in the final results of our 
separation tests. Hence, the presence of absorptions does not 
depend on our analysis of the continuum or on the editing 
'style' of the mean spectrum during the separation, nor does 
it depend on the choice of iterative sequence and the system of 
RVs we adopt while performing the separation. 

Despite what is stated above, and in order to strengthen 
our confidence in the result of the separation, we measured 
the equivalent widths (EWs) of the 5 absorption lines (i.e. in- 
cluding this time He) seen in the O spectrum. If the separation 
process is reliable, we might expect that these EWs should not 
depend on the orbital phase. We also measured the EW of the 4 
most prominent emission lines in the WR spectrum. The results 
are shown in Fig. [5] It appears that the EWs of the absorption 
lines in the individual extracted O-star spectra are roughly con- 
stant. No obvious periodic pattern or phase-locked variations 
are detectable in these EWs. Moreover, the EWs of the emis- 
sion lines in the WR spectra are also relatively constant (see 
Fig-H]- bottom figure). 

Continuing to test the sensitivity of the overall appearance 
of the weak absorption components to details of the extraction 
algorithm, we attempted to construct an unperturbed WR spec- 
trum. For that, we shifted the WR spectra to the same reference 
frame using the RVs of the accepted orbital solution. Then, we 
produced an average from the spectra around phases 0.4-0.6 
(corresponding to inferior conjunction of the O star). The red 
side (up to 2000 km s ) of the profiles was assumed to repre- 
sent an unperturbed profile. Then, we combined the 3 spectra 
around phase 0.15 and considered the blue side of the average 
(down to -2000 km s _1 ) as a fair representation of the unper- 
turbed WR line profile. Then we merged the two sides for each 
line and subtracted the templates from the individual spectra. 
Unfortunately, the relatively low S/N ratio of the spectra dis- 
allowed any further analysis of the absorption features, which 
were submerged in noise in all individual spectra. 

A montage of the most prominent lines in the WR rest 
frame can also be used to understand more precisely what 
happens. The plots of H<5, Hell /14200, Hy (all coming from 
the EMMI data) and Hell /14686 (our dataset-1) are shown in 
Fig. [6] The montage of the Hell /14686 profiles shows very 
clearly the WWC excess-emission, traveling from the red side 
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Fig. 5. Top figure: EW of absorption lines measured on the ex- 
tracted O-star spectra. Bottom figure: EW of emission lines 
measured on the extracted WR spectra, including the WR ab- 
sorption features. Note that panels of the two figures do not 
have the same vertical scale. 

(around phase zero) when the O star is behind and the shock 
cone is pointing away from the observer, and reaching its blue- 
most position around phase 0.5. On the montages of the three 
other lines, it can be seen that the absorption feature is most 
clearly seen between phases 0.15 and 0.25. This is consistent 
with the fact that an excess emission is present in the red-center 
portion of the profile during that narrow phase interval. Then, 
when traveling to the blue, around phase 0.5, this emission ex- 
cess fills in the absorption feature, then retreating to the red side 
at phases around 0.9 and producing a profile with a pronounced 
redshifted bulge. As mentioned before, the latter caused prob- 
lems during the separation process. 

The ultimate test to check if the WR absorption lines do be- 
long to the WR star is to measure the blueshifts of these absorp- 
tions, assessing both the blueshift on the mean WR spectrum 
(that is a measure of the separation between the WR absorp- 
tions and emission lines) and also blueshifts seen in individual 



lines (to prove that the blue-shifted absorption lines follow the 
WR orbital motion). To measure the blueshifts on individual 
spectra, we first shifted the WR spectra to the same frame us- 
ing the orbital RVs, then combined them in orbital phase bins of 
2 or 3 spectra in order to increase the S/N ratio. A double gaus- 
sian fit was performed on the 4 prominent lines cited above: 
Hell ^4025, H<5, Hell A4200 and Hy, hence providing a posi- 
tion both for the absorption and the emission profile. The re- 
sults are displayed in Fig. (including only the measurements 
for which the fit provided a reasonable solution). The squares 
show the velocities of the absorption lines compared to the rest 
wavelength. The open circles show the velocities of the emis- 
sion lines. 3 




Fig. 7. Velocities of the WR absorption and emission lines ob- 
tained from the positions of the double gaussian fits, versus or- 
bital phase. The positions are compared to the rest wavelength 
of the corresponding line. For the H+He blends, the simple av- 
erage of the rest wavelengths has been taken. Filled squares: 
absorption lines. Open circles: emission lines. The crosses in- 
dicate the orbital RVs of the O-star in the WR frame. 

Velocities of absorption and emission lines in Fig. Q show 
a rough trend to be slightly in anti-phase. This is particularly 
clear for Hy where initially the RVs of emission and absorption 
are well separated, then become almost identical at phase 0.5, 
and again move apart at phase 0.9. Recall that all measurements 
are made in the WR rest-frame. Clearly, the excess emission 
coming from the WWC dominates the fits, apparently shift- 
ing the absorption line. However, the amplitude of the shifts of 
the absorption lines is much smaller than the corresponding O- 



3 For the H+He blends, the simple average of the rest wavelength 
has been taken. The important point here is to compare the veloci- 
ties of the absorption and emission profiles together (see also previous 
footnote). 
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Fig. 6. From left to right: Montage of H6, HeII4200, Hy and HeII4686 from the WR component plus the WWC emission. All 
spectra were shifted to the WR rest-frame. For HeII4686 we used a reduction scale factor of 15. In the first 3 figures, the mean 
position of the absorption line is indicated with a dashed gray line. 



star RVs, whose amplitude would be "enlarged" in the WR rest 
frame (the O-star motion is indicated by crosses in Fig.0, and 
thus the absorption lines cannot belong to the O component. 
Disregarding the WWC-induced drift, the blueshift of the ab- 
sorption lines can therefore be considered as constant relative 
to the rest wavelength, in the WR frame. 

Finally, we took all the shifted WR spectra except those at 
phase 0.9 and combined them into a simple average, to increase 
the S/N ratio. We then performed a double-gaussian fit on the 4 
prominent lines. The results are summarized in Table |4] where 
the positions of the absorption and emission lines are compared 
to each other. The blueshift is clearly increasing with lower 
B aimer/Pickering number(s). 



Table 4. Blueshifts of WR absorption lines when compared to 
the positions of the emission lines. 



Line(s) 


Blueshift (km s~') 


Hel ,14025 + Hell ,14026 


-86.3 + 18.3 


m + Hell ,14 100 


-161.8 + 12.5 


Hell ,14200 


-228.7 + 11.5 


Hy + Hell ,14338 


-289.0 + 10.6 



Altogether, these tests on the absorption lines seen in the 
WR spectrum allow us to claim that these blue-shifted lines be- 
long to the WR star itself. Therefore, they are probably formed 
in the wind of the WR star, similarly t o what has been seen in 
other single SMC and LMC stars (see lFoellmi et all2003al in 
particular). 



4.2. Spectral types of the WR and O stars 

Although they have a lower quality and the dataset is not ho- 
mogeneous, we attempted to separate the orbital components 
with the spectra of dataset 1, using a blind separation with 
5 iterations, the orbital RVs and no editing. Thanks to the 
broader wavelength coverage and the lower resolution, we ob- 
tained a reasonably good O-star spectrum. The result is dis- 
played in Fig. [8] A comparison of the mean O-star spec- 
trum in BAT9 9 129 with standa r d spe ctra in the same (blue) 
range of IWalborn & Fitzpatrickl £l990) leads to a type 05V. 
More precisely, the most frequently used features for deter- 
mining the spectral class are Hell ,14541 and Hel /14471. In 
our spectrum, the former is clearly stronger than the latter, 
pointing to a spectral type earlier than 07 (see close-up in 
Fig. [S}. This is confirmed by the comparison of Hell /14200 
and Hell + Hel ,14025. We also measured the equivalent width 
of the lines Hell /14541 and Hel ,14471 to compute their ra- 
tio. F ollowing the spectral classification relation from Mathvs 
dl988l) . we obtain log[£W 44 7i/£W4542] = -0.48,-0.42 de- 
pending on the continuum position, with a preference for the 
first value. This corresponds to spectral types 05 and 05.5 re- 
spectively. Nonetheless, given the noisy nature of the spectrum, 
even a spectral type 06 cannot be excluded. As for the lumi- 
nosity class, the extracted absorption lines are relatively broad, 
with no trace of SilV beside H6 on the higher-resolution EMMI 
spectra, hence broadly suggestive of a V-III luminosity class. 
Moreover, the line Hell ,14686 in Fig.[S]is (particularly) strong 
in absorption, and this defines a luminosity class V. 

The WN spectrum resulting fro m the separation loop 
can b e compared to the criteria of Smit h. Shara. & Moffat 
(1996). Using the peak values, we obtain a ratio for 
Hell ,15411/HeI ,15876 = 16.1. This points towards a WN3 
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Fig. 8. O-star spectrum obtained by a separation process on the 
spectra of dataset 1 . The inset shows the region around the two 
lines Hell AAA! 1 and Hel /I4542 used for classification. 



subclass, that is confirmed by the peak ratios of NV /14604 
and NIII /14640, since the latter is very weak. The practical ab- 
sence of the CIV /15808 line makes the result definitive: WN3. 
Finally, we measured the Pickering emission -line decrement. 
Using the description of ISmith et alJ dl996l) . the heights of 
H+He lines clearly exceed the line drawn between pure Helium 
lines. Results are shown in Fig. [9] Moreover, the two hydrogen 
criteria of Smith et al. are equal to 0.22 and 0.24, respectively. 
These are probably lower limits since the H-containing blends 
are more suppressed by stronger absorption components than 
the pure He lines. Combined with the presence of H-absorption 
lines, this leads to the moderate presence of hydrogen and a 
final type of WN3(h)a. Therefore, BAT129 appears to be an 
early-type hydrogen-containing WN star with a dwarf 05 com- 
panion. 



4.3. Luminosity ratio 



Using the methods described in ICher enashch uk et"afl Jl995h 
we derived the(visual) luminosity ratio q = Lwr/Lo- To 
achieve a meaningful ratio, we combined EMMI spectra be- 
tween phases 0. 1 1 and 0.25, i.e. outside the eclipses phases (the 
spectra of dataset 1 were not of sufficien t quali ty and resolution 
for this). The classical method by Beals ( 1944) provides an up- 
per limit, and the new method by Cherepa shchuk et all (^995) 
gives the "correct" value (i.e. corrected by the fact that the ab- 
sorption and emission lines are superimposed on each other). 
The determination of the luminosity ratio is made by com- 
paring the equivalent widths (EW) of O-star absorption lines 
both in the original mean WR+O spectrum and the "pure" O- 
star mean spectrum. We measured the EWs on the following 
lines, by gaussian fitting: He, Hell A4025, R6 + Hell ,14100, 
Hell ,14200 and Hy + Hell ,14338. 




0.9 



4200 4400 4600 4800 5000 5200 5400 
Wavelength [A] 

Fig. 9. WR spectrum obtained by a separation process on the 
spectra of dataset 1. The line is used to show the Pickering 
emission-line decrement, leading to a moderate (h) presence of 
hydrogen in the WN spectral type. 



Two problems arise here. First, the absorption lines belong- 
ing to the WR star in the WR+O spectrum must be removed. 
To account for this, we corrected the EW of absorption lines in 
the original WR+O spectra by measuring the EW of the WR 
absorption lines in the extracted WR spectra: 



EWbin = EWwR+O - EW a b s WR- 



Also: 



EW. 



single 



EW bi n ■ LCR 



- 1, 



(1) 



(2) 



where LCR stands for "Line-to-continuum ratio" and is the 
measure of the line strength to th e continuum level (i.e. peak 
value) in the WR+O spectrum (seelcherenashc huk et alJl995l 
for the detailed equations). 

The second problem is to compare our values of EW S i„ g i e 
to that of single O stars. We first measured the EW on the 
extracted O-st ar spectrum. Then, c omparing our results with 
the values of IConti & Frosil i 19771) for the lines Hell /14025, 
Hell /14200 and Hy, it seems that our values for the helium 
lines are overestimated, while that of Hy is underestimated (by 
about 0.3 and 0.6 A respectively). Taking lower values for the 
Helium lines give inconsistent results, since the EW on the 
combined binary spectrum would be larger than that of the sin- 
gle spectrum. Taking the upper value for Hy would actually 
correct the result of this specific line to a level comparable to 
that of the other lines. Although not fully reliable (since the re- 
sult depends heavily on the output of a complicated separation 
process), we consider that our measurements on the extracted 
spectra are nevertheless globally consistent. 
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Assuming of course that there is no third-light contamina- 
tion, and taking these values with care, we obtain the following 
mean luminosity ratios: 0.39 + 0.18 and 0.28 + 0.20, respec- 
tively with the two methods. In the following, we assume a 
luminosity ratio of 0.3. 

4.4. The O-star characteristics 

Given the spectral type we assigned to the O star, the 
temp erature-spectral type relation published by iMassev et alJ 
(2005) in their detailed study of LMC O stars, implies an effec- 
tive temperature of about 40 kK, interpolating roughly the val- 
ues between the SMC and the Milky Way relations. Similarly, 
using t he scale for a solar m etallicity and for galactic dwarf 
stars by Marti ns et alJ d2005l) implies T e ff = 41 kK for a spec- 
tral type 05V. However, these latter authors emphasize metal- 
licity effects which provide higher temperature at lower metal- 
licity. If we consider a shift of about one spectral subclass, this 
gives (for an 04V star): Teff = 43 kK. For such a temperature, 
the luminosity of the O star is log(L/L ) = 5.49. Let's assume 
this value of L for the secondary star in BAT 129. 

One of the stars in the sample of IMassev et alJ J2005), 
namely LH 58-496, has the same spectral type (05V) and thus 
the assumed luminosity (log(L/L Q ) = 5.49) of the O star in 
BAT129. The radius of the star quoted by the authors is 10.5 
R . Following EggletonJ dl983l) . this value is slightly smaller 
than that of the Roche radius of the O star in BAT 129 (11.8 
R Q ). Moreover, the evolutionary mass of LH 58-496 is 40 
M , compared to the measured Keplerian mass of BAT129: 
M sin 3 i = 23M (see Table QJ. To reach a value of 40 M 
for BAT129, the inclination angle of BAT129 should be 57 de- 
grees. Given the strong eclipsing nature of BAT129, this value 
of the inclination should be considered as a minimum estimate. 
Considering the large radius of the 05 star, the secondary com- 
ponent of BAT129 well may be an 06, with a smaller radius. 

In their sample, Mass ev et al ] J2()()5l) have an 06.5V star, 
AV 446 (in the SMC though). It has a temperature of T eff = 
41.0 kK, an evolutionary mass of 33 M , and a radius of 8.8 
R (i.e. 75% of the O-star Roche radius). These values are also 
quite consistent with the estimates we draw on the secondary 
star ofBAT129. 

Given the short period, it is likely that the system is syn- 
chronized. If so, then with a radius of 8.8 R Q , the equatorial 
velocity becomes 161 km s . However, the FWHM measured 
on the extracted O-star spectrum gives a projected rotational 
velocity slightly larger by about 20 to 40 km , as shown in 
Table |3 (assuming that the dominant broadening cause is the 
Doppler effect). On the other hand, these values are very close 
to the synchronized velocity of the star if it had a radius of 10.5 
R Q (192 km Oversynchronicity is alread y known to occur 
in so me WR+O bin aries, such as V444 Cy g ( March enko et alJ 
1 1994 and CQ Cep jMarchenko et alJl995i) . 

4.5. The WR-star characteristics 

The visual app a rent magnitude of BAT 129 quoted by 
iBrevsacher et alJ i 19991) is 14.92, refering to a work by 



Table 5. Projected rotational velocities inferred from the O- 
star absorption lines, assuming that the broadening of the lines 
is due only to Doppler broadening. Values have been corrected 
by the 2.6-pixel resolution. 



Line 


FWHM (A) Velocity (km s" 1 ) 


He + Hell ,13968 


6.72 


195 


HeH ,14025 + Hel ,14026 


5.51 


158 


H<5 + Hell ,14100 


7.43 


209 


Hell -14200 


6.61 


181 


Hy + Hell ,14338 


7.19 


191 



Rousseau et al. (1978). Assuming a distance modulus of 18.5 
for the LMC ( Ivan d en Bergh 2000) and IS visual extinction of 
Av = 0.4 (i.e. 3.1 * 0.13 see lMassevll2003l IConti & Morril 
1990), the absolute visual magnitude of the system is -3.98, i.e. 
fainter than the magnitude o f the sing le 06.5V star mentioned 
above (Mv = -4.7, Massev et al. 2005). The Johnso n magnitude 
ofBAT129 given bv lFeitzinger & IsserstedlH 1983b is 14.68, i.e 
Mv = -4.2, which is still fainter than the O star value 4 . Two ob- 
vious reasons can explain the above discrepancy. First, the mea- 
sured magnitude of the system is too faint because the eclipses 
were not accounted for. This effect is unlikely to resolve the 
discrepancy, since it is smaller than 0.2 magnitude (see Fig. [2}- 
Second, the spectral type inferred from our extraction process 
is wrong and a cooler spectral type must actually be assigned 
to the O star in BAT129. This second possibility is more plau- 
sible, although we have no definite means to verify it. 

We can nonetheless compare the properties of the 
WR star in BAT 129 with LMC WN stars studied by 
lHamann & Koesterkd J2000I) . Taking our luminosity ratio of 
0.3, the WR star has a luminosity of log(L/L ) = 4.97, and 
M^-3.34. This absolute magnitude is nonethel ess close to th e 
average value for the Galactic WN3 stars JvanderHuchtl200ll) . 
and slightly fainter than that of WR3 in our Galaxy (Mv = -3.7, 
iMarchenko et all2004l) 

Given the luminosity of the WR star in BAT129 and its ab- 
solute visual magnitude , the most similar star in the selected 
subsample of lHamann & Koesterkd l l2000l) is BAT99 63 (Brey 
52): it has a spectral type WN4, a temperature of 71 kK, 50% 
hydrogen (by mass), a radius R* (i.e. at a Rosseland optical 
depth of 20) of 4.7 R and a mass-loss rate log(M/M ) = 
-5.2. This mass-loss rate value is i n agreement with the re - 
cent theoretical value computed by IVink & de Koterl J2005), 
for the metallicity of the LMC: -5.26, -5.49 for /3 (exponent 
of the wind velocity law) equal to 1 and 3, respectively. In 
the MACHO light curve, spanning more than 10 years, we 
have also searched for any significant period changes. We have 
split the lightcurve into 10 equally long bins, and performed a 
Fourier analysis on each of them. Note that given the number 
of datapoints and their spread in time, the maximum number of 
bins allowing a precise analysis is not much larger than 10. We 
found no significant period changes, with a precision of 0.003 
days over the whole timespan. However, this precision does 

4 We were not able to compute the correct MACHO magnitude, fol- 
lowing equations 1 and 2 of lAlcock et al. ( 1999), since some coeffi- 
cient values were not accessible on the Internet, at the time of writing. 
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not allow us to put any meaningful constraint on the WR mass- 
loss rate. Comparing the Keplerian mass of the WR component 
in BAT12 9 (M^sin 3 i = 14 M ) t o the actual mass of BAT63 
quoted by Hamann "& Koesterkd d2000l) - 22 M - the inclina- 
tion angle in BAT 129 must be 59 degrees, which is equivalent 
to the estimate drawn from the O star. 

For the moment, let us assume that the WR star in BAT129 
has followed a pure single-star evolution. The WR star must 
have had an initial mass above that of the current secondary 
star, i.e., at least 35 M . On the other hand, if we assume that 
the WR star is, say, in the middle of its eWNL phase (i.e. 
its h-rich WN phase, see below), we can use the mass loss 
rate quoted above to estimate the amount of mass lost by the 
star during the WR stage. The duration of the Helium-burning 
phase of a 40 M star, followin g the model s of rotating WR 
stars at Z=0.008 of lMevnet & Maederl J2005h . is 0.550 million 
years. Assuming the WR star in BAT129 is well described by 
these models, this means that the star has lost slightly less than 
2 M during the WR phase. The 60 M model would imply 
an even lower value, but even the 30 M model does not give 
a value larger than 3 M lost by the WR. Taking the slightly 
smalle r mass-loss rate of W R3 in our Galaxy (log(M/M ) = 
-5.62, iMarchenko etail 2004) of course does not resolve the 
situation. This means that the mass-loss rate quoted above is 
certainly not valid for the WR star in BAT129 as a single star, 
since the estimated mass lost is too small by a factor 5 (2-3 M 
compared to at least 10 M from its initial mass). 

With the assumption that BAT 129 is a non-interacting sys- 
tem and that no RLOF occurred in its past evolution, some of 
the parameters of the WR star and its companion can be rel- 
atively well explained by single-star evolution. However, the 
combination of the possible oversynchronous rotation of the 
O star, the very short period of the system and the theoretical 
single-star mass loss rate clearly too small seem to point toward 
an interacting scenario for BAT 129, although the presence of 
hydrogen in the WR star might be contradictory at first sight. 
This is discussed in the next section in a broader perspective of 
the WN evolution with hydrogen. 

5. Discussion 

It appears that weak absorption features are superposed on the 
emission line profiles in the extracted WR spectra. We have 
shown that the features are moving with the WR star. These ab- 
sorption features are blueshifted relative the main emission pro- 
files, and the blue shift seems to increase in the lines with pro- 
gressively lower Pickering (or Balmer) number. Therefore, we 
are tempted to conclude that these absorption features are iden- 
tical to the intrinsic absorption lines seen in many early-type 
WN Magellanic Cloud stars and WR3 in our Galaxy. However, 
for the first time we detect them in a short-period binary. The 
presence of hydrogen in an early-type WN star in a binary has 
important consequences that are discussed below. 

5.1. Wolf-Rayet star classification 

The presence of hydrogen in WN stars is relevant for their clas- 
sification and from the point of view of their evolutionary sta- 



tus. Most of the late-type WN stars (WNL) have hydrogen in 
their winds, as they are believed to be the first WR stage after 
the main-sequence and/or a possible transition stage. However, 
a significant number of exceptions exist to this rule. Similarly, 
most of the early-type WN stars (WNE) are hydrogen-poor, or 
hydrogen-free, but again, with a significant number of excep- 
tions. For some reasons, that are likely to be due to rotation, 
the exceptions are more numerous in the SMC. Five out of 1 1 
WN stars in the SMC have a very hot spectral type (WN4, 
WN3) and present signatures of a large amount of hydrogen 
(ratio H/He by number ~ 1, P. Crowther, private communi- 
cation). Therefore, the ionization subclass (i.e. the tempera- 
ture) is not the decisive criterion to distinguish between late 
and ear ly-type WN stars. This is the reason whv lFoellmi et alj 
(2003b) proposed a complementary and evolutionary classifi- 
cation for WN stars (recognizable by putting an "e" in front of 
the class) which is based on the hydrogen content. The main 
eWNL/eWNE distinction is now based on the amount of hy- 
drogen. To discriminate the h-rich from the h-free stars, the 
"b" (for "broad" lines) classificatio n criterion of the classifi- 
cation scheme of Smith et alj II 19961) can be used. As a matter 
of fact, Hama mTet al.Nl995l) have shown that the width of the 
emission line s is directly conne cted to the hydrogen content 
(see also ISmith & Maeder ll998h . Moreover, since a clear anti- 
correlation between the "b" and "h /(h)" spectral not ations of 
WN stars is observed in L MC stars dFoellmi et al.l200 3b) - an 
effect already mentioned by Smit h et alT7l996l) - the "b" crite- 
rion is an "antipode" of the "h" criterion. 

This classification was meant for single stars, since the hy- 
drogen content in WN stars in binaries was not known. With 
BAT129, we can extend the applicability of this classification. 
It is rather natural to think that the same causes produce the 
same effects: when the lines are narrow, it means that hydrogen 
is present in the atmosphere of a WR star, and in some cases 
where the amount of H is substantial, an alterning Pickering 
emission-line serie and/or absorption lines may be visible in 
the optical spectrum. Moreover, we note that this classification 
implies that the SMC contains virtually no eWNE stars, but 
only eWNL stars (and one WC/WO binary). This fact is well 
reproduced b y the recent model s of ro tating WR stars at low 
metallicity of M evnet & M aeder (2005, see their Fig. 9). 

5.2. Evolution of WR+O binaries 

A simple estimate of the minimum period of BAT129 allows us 
to show that the system was certainly interacting in the past. Let 
us assume a ZAMS mass of the O star similar to the estimated 
mass above (say, 35 M Q , i.e. the star is assumed to have lost 
no mass), and t he corres ponding Main Sequence (MS) radius: 
~10R Q (see e.g. lMartins et all 2005. their Table 1). Assuming 
that the present WR star had an initial mass of, say, 45 M , 
and became a WR star by simple wind mass-loss, its maxi- 
mum radius during the MS phase is the radius of a 45 M MS 
supergiant: ~20R Q (i.e. not conside ring a red supergian t nor a 
Luminous Blue Variable phase; see Mart ins etalJ koOS. their 
Table 3 for the radius value). Assuming that the minimum dis- 
tance between the two stars is at least the sum of their radii 
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(dmin > 3O/? ), and using Kepler's third law, the minimum pe- 
riod of the system is 2.2 days, which is very similar to the cur- 
rent period. Therefore, given the extreme example considered 
here, the widening of the orbit from 2.2 to 2.7 days caused by 
a simple mass-loss by the stellar wind during the WR phase 
of the primary component can be ruled out. This implies that 
to obtain the present short period, the system must have gone 
through an interacting phase that shrank the orbit. 

Furthermore, it is very interesting to compare the actual pa- 
rameters of BAT 12 9 with the latest models of binary star evolu- 
tion of Pe trovic. Langer. & van de r Hucht ( 2005a). The models 
of binary stars that lead to WR+O systems show that Roche- 
lobe Overflow (RLOF) and mass transfer occur for short- 
period systems. We can distinguish Case A mass transfer, that 
starts during the core hydrogen-burning phase of the primary, 
and Case B mass transfer starting during the shell hydro gen- 
burning phase. As emphasized by Pe trovic et alJ (|2p05a), the 
net effect of a Case A (that also contains an intermediate Case 
AB mass transfer episode) or a Case B scenario is the in- 
crease of the period. If such mass transfer occurred in the 
past in BAT 129, this means that the initial period would have 
been even smaller than about 2 days, which is very unlikely. 
Therefore, following Petrovic and collaborators, the evolution 
of BAT 129 should be considered as the probable result of con- 
tact evolution. 

IWellstein. Langer. & Braunl J200ll) have studied the for- 
mation of massive contact binaries. Unfortunately, among the 
74 models they have computed, none of their models actu- 
ally include primary stars more massive than 25 M Q . If we 
nonetheless consider globally their results, their models also 
lead to an increase of the period after the RLOF (whatever its 
kind). We note here that an oversynchronous rotation veloc- 
ity of the secondary star is known to occur in other WR+O 
syste ms, among which is th e famous Galactic syste m V444 
Cyg llMarchenko et alJll994l) . WR127 lMassevHl98ll) and CQ 
Cep ( Marchenko et al. 1995). Interestingly, while CQ Cep is 
supposed to have experienced a contact phase, this is not the 
case for V444 Cyg, which is relatively well explained by a sta - 
ble mass transfer scenario, according to lPetrovic et alJ J20 05a). 
The uncertainty on the radius of the WR star in BAT 129 is un- 
fortunately too large to allow a conclusive discussion on this 
phenomenon related to the interaction of the two stars. 

Altogether, it seems likely that BAT 129 followed an inter- 
acting (possibly contact) evolution. A detailed analysis of the 
abundances of the O star could possibly help to probe the ques- 
tion of whether the companion star has accreted mass from 
the primary, and how much. We have obtained very high S/N 
FORS 1 spectra at the VLT of BAT 1 29 that could possibly allow 
such analysis. This is left for a forthcoming paper. 

5.3. Gamma-ray burst progenitors 

In this section we try to consider the binary star BAT 129 or 
its WR component alone as "test particles" for the scenario of 
formation of Gamma-ray bursters (GRB). 

Although WR stars are believed to be the progenitors 
of Type Ib/c supernovae associated with long-duration GRBs 



quite successfully described by the "col- 
lapsar" model (MacFa dven & WooslevlH99^) . it is still very 
difficult to answer the questi on: what kind of WR stars are the 
best can didates for a GRB? IPetrovic. Langer. Yoon. & Heg erl 
J2005bl) try to address this issue from a theoretical point 
of view. They emphasize that three essential ingredients are 
needed to produce a GRB: a massive core to form a black hole 
(BH), sufficient angular momentum to form an accretion disc 
around the BH with bipolar jets, and loss of hydrogen envelope 
to let bipolar jets reach the star surface. They show that both 
single and binary non-magnetic WR stars can retain enough 
angular momentum for a GRB to form. Hence the question re- 
mains: are binaries more favorable GRB progenitors than sin- 
gle stars? 

The recent mo dels of single rotating WR stars at differ- 
ent metallicities of lMevneT& Maederj J2005t) show 3 impor- 
tant things. (1) Final surface velocity at low metallicity is 
higher than at high metallicity, a fact that is explained as a 
consequence of the lower mass loss rates at lower metallic- 
ity. Therefore, low metallicity galaxies such as the SMC and 
the LMC are better locations to find fast rotating WR stars 
at the en d of the He burning p hase, in agreement with the 
results of iLe Floc'h et al ] fc003h who show that GRB occur 
mostl y in blue sub-luminous galaxies (see also lProchaska et all 
2004). (2) At a given metallicity the velocities obtained at 
the end of the He-burning phase are larger for smaller initial 
mass stars. (3) The last point is that most of the models of 
iMevnet & MaedeJ J2005b at low metallicity have enough spe- 
cific angular momentum in their core to be good candidates 
for collapsar models. Again we can ask: are binary stars bet- 
ter GRB candidates than single stars, even if single star WR 
evolution provides the necessary conditions for such an event? 

How does BAT 129 fit in this picture? While the metal- 
licity is that of the LMC, and the combination of the very 
early spectral type and t he presence of hydrogen is probabl y 
a sign of rapid rotation JFoellmi et aljboo"^ iFoellmll l2004) . 
the main problem is to know what was really the initial mass 
of the WR star. If it had followed a single-star evolution, it is 
likely that it was initially a lower-mass star (say with an ini- 
tial mass between 20 and 40 M Q ). This is consistent with its 
la rge rotational velocity (s ee again the corresponding model 
of Mevnet & Maeden l2005). and the clear presence of hydro- 
gen. However, we have shown above that the characteristics of 
BAT 129 as a binary make it a good candidate for an interact- 
ing binary, possibly with a contact phase. We can therefore ask: 
could the initial mass of the WR star be very large (say above 
50 M ), and most of it has been lost in the interaction of the 
binary? Could binary evolution produce such a WR+O system 
where the WR star looks like its SMC counterparts? If the WR 
star lost a large amount of mass and the secondary did not ac- 
crete all the material, it is likely that a nebula should be seen 
around BAT129. We find however no clear nebular lines in our 
spe ctra. 

Tutukov & Cherepashchuk (2004) also mention that the 
creation of a BH and the generation of a GRB associated with it 
is most probable in the closest massive binary systems (not only 
WR+O but also WR+ [compact companion] binaries), since the 
presence of a close companion limits the rotational angular ve- 
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locity of the WR core to the orbital angular velocity of the bi- 
nary system. Very short period WR+O binaries formed through 
a common envelope phase, as BAT 129 could b e, is very rele- 
vant in this picture (see also the models bv lWellstein & LangeJ 
Il999l) . Although not a definitive argument, we note that the rar- 
ity of very close WR binary systems might help the frequencies 
of GRBs and fast r otating WR core collapses to agree (see e.g. 
IPostnov & Cherepashchukl200ll) . 

The point of GRB progenitors is further investigated by 
Hirschi. Mevnet. & Maede for single star evolution. 

They note in particular that angular momentum at the end of 
the He-burning is a good approximation of the angular momen- 
tum contained in the core at the collapse. This is an important 
point for observational tests, since even the pre-core collapse 
evolution phases are rather short and unlikely to be extensively 
observed. The authors also give a rough estimate of the initial 
mass range to form a BH, that is between 20 and 60 M . But 
forming a BH is not sufficient to produce a GRB; the angular 
momentum must be large enough to build up an accretion disk 
that will power bipolar jets. Interestingly, they show that WR 
stars of the WO subtype are the most probable GRB progeni- 
tors, since they are for med at low metallicity, and always form 
a BH (see also lPostnov & Ch erepashchuk 2001). However, the 
initial mass of a WO star is estimated to be > 5OM . 

All the above considerations lead to the reasonable conclu- 
sion that GRB progenitors are certainly a special trade-off of 
the different parameters of WR evolution. Some things seem 
however to be clear. The metallicity should be low, and the 
initial mass of the WR star not too high to avoid any strong 
mass loss, and in particula r any Lum inous Blue Variable (LBV) 
event (see IMevnet & Maede rl2005l for limiting masses of the 
3 main evolution sce narii of W R stars). I t is int eresting to note 
that in the models of Mevnet & MaedeJ (|2005) at SMC metal- 
licity (Z = 0.004), only the 60 M model has clearly a higher 
rotational velocity at the end of He burning, while the low- 
est initial mass models (30 and 40 M ) at LMC metallicity (Z 
= 0.008) have the largest rotational velocity. We may wonder 
whether only single WO stars at the lowest metallicity produce 
GRBs, while lower mass WR stars do so at slightly higher 
metallicity. In that context, the unique WO binary star of the 
SMC (WR8) deserves probably much more attention (see e.g. 
iBartzakos et all200lh . 

One may also want to discuss the growing amount of ev- 
idence that hydrogen is clearly detected in GRB afterglows. 
This has been first discussed by Rigon et al. (2003) for super- 
nova SN 1999E possibl y associated with the BATSE event of 
GRB 980910. Recently. Istarling et alJ J2005t) found the clear 
signature of hydrogen in the spectrum of t he ear ly afterglow 
of GRB 021004. Finally, Ivan Marie etall J2005I) have made 
evolutionary models of circumstellar material around a mas- 
sive star that produces a GRB. It is therefore relevant to search 
for the characteristics of hydrogen in WN stars. As suggeste d 
by the analysis of WR3 in our Galaxy (Mar chenko et al. 2004), 
the absorption lines form relatively close to the stellar core. It 
is likely that the absorbing material is then gradually acceler- 
ated to the terminal wind velocity. Since we have shown that it 
is not unusual to see a significant amount of hydrogen in single 
and binary WN stars at low metallicity, it might be not too sur- 



prising to observe high-velocity hydrogen absorption f eatures 
in GRB afterglows, such as described by Starli ng et alJ J2005) 
on GRB 021004. 

On the other hand, van Ma rie et all J2005I) argue that the 
presence of an intermediate velocity component in the after- 
glow of GRB 021004 implies that the WR phase was short, i.e. 
that the WR shell was still intact when the star exploded. This 
implies that the initial mass of the progenitor must have been 
small (i.e. about 25 M ), since the s maller t he initial mas s, th e 
shorter the WR phase (see Fig. 9 in IMevnet & Maederll2005l) . 
But such stars have a very short H-free WC stage. This is pos- 
sibly in contradiction with th e requirement of no hydrogen left 
in the atmosphere (MacFadven & Wooslev 1999). 

It is interesting to note that theoretical models seem to pro- 
duce two very different t ypes of GRB progenitors. On one 
hand, iHirschi et alJ d2005l) show that single WO type stars are 
the best GRB candidates, while the presence of hydrog en in 
GRB afterglows and the results o flvan Marie et alj (|2005) point 
toward lower mass WR stars. Could these latter stars still pro- 
duce GRB if they belong to an interacting binary system, where 
the companion helps to remove the hydrogen layers? The ques- 
tion remains open. 

6. Conclusion 

We have presented new medium-resolution NTT-EMMI spec- 
tra of the LMC eclipsing Wolf-rayet binary BAT99-129. We 
have been able to separate the spectra and find a reliable or- 
bital solution. We assign a spectral type WN3ha + 05V to the 
system. The O star seems to appear possibly oversynchronous 
with the orbital motion. We find that the WR star contains a 
detectable amount of hydrogen, seen in absorption in its spec- 
trum. While its actual characteristics are shared w ith the eWNL 
stars seen in the SMC ( Foell mi et alj2003atlFoellml2004i) . the 
evolution that led to the present state is certainly much differ- 
ent. Interacting and binary models that could be applied to sys- 
tems like BAT 129 (hence also to V444 Cyg) would be of major 
interest. 

BAT99-129 appears to be a very interesting piece of the 
puzzle of GRB progenitors. However, the most important ques- 
tions remain: What are the kinds of the progenitors of the long 
soft GRB progenitors? Are they single or binary stars? Are they 
low mass or high mass WR stars? The presence of absorption 
lines in BAT 129 proved to be of the greatest interest and fur- 
ther studies of GRB afterglows will certainly confirm the im- 
portance of these absorptions. 

BAT 129 is definitely a very interesting target for X-ray 
observations, considering the eclipsing nature of the system 
and the clear signs of wind-wind collision. Unfortunately, we 
found no XMM, ROS AT, nor Chandra observations of this sys- 
tem. Moreover, the eclipses deserve much more study. In par- 
ticular, the eclipses can be used to probe the actual sizes, in 
physical units, of the line emitting regions. This is now possi- 
ble thanks to the complete orbital solution. We obtained VLT 
FORS 1 spectra for this purpose and their analysis will lead to 
a forthcoming paper. 

All in all, BAT99-129 appears to be very similar to the well 
known WR binary in our Galaxy, V444 Cyg. 
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